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Abstract 

In this paper we study scalar multivariate subdivision schemes with general 
integer expanding dilation matrix. Our main result yields simple algebraic 
conditions on the symbols of such schemes that characterize their polynomial 
reproduction, i.e. their capability to generate exactly the same polynomials 
from which the initial data is sampled. These algebraic conditions also allow 
us to determine the approximation order of the associated refmable functions 
and to choose the "correct" parametrization, i.e. the grid points to which the 
newly computed values are attached at each subdivision iteration. We use 
this special choice of the parametrization to increase the degree of polynomial 
reproduction of known subdivision schemes and to construct new schemes 
with given degree of polynomial reproduction. 

Keywords: multivariate scalar subdivision schemes, expanding dilation ma- 
trix, polynomial reproduction, refinement parametrization 

1. Introduction 



In [2j, the authors derive simple algebraic conditions on the subdivi- 
sion symbol of a multivariate scalar subdivision scheme with dilation ma- 
trix M = ml, m G Z, \m\ > 1, that allow us to determine the degree of 
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polynomial reproduction of such a scheme and, consequently by 16|, its ap- 
proximation order. In this paper we extend the results of [2J to the case 
of general expanding dilation matrix M £ Z sxs whose spectral radius p(M) 
satisfies p(M) > 1. On the one hand, our interest in the study of such a gen- 
eral dilation matrix M is motivated by the existence of multivariate scalar 
subdivision schemes for which the results of ^ are not applicable, e.g. 

M = ( n J , a, b, c £ Z, \a\, |c| > 1, 
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101 . |14| . On the other hand, even for schemes whose dilation matrix 
M satisfies M e = ml for some integer i > 1 and m £ Z, \m\ > 1, e.g., 
y/3— subdivision with M 2 = —31, the results of this paper allow us to di- 
rectly investigate the properties of the associated symbol, instead of dealing 
with its £— th iterated version. 

Moreover, in contrast to polynomial generation, i.e. the capability of a 
scheme to generate the full space of polynomials of certain degree, there 
are very few theoretical results on polynomial reproduction of subdivision 
schemes, (see |, B El). Thus, we strongly believe that our paper may 
further contribute to the development of this topic. See ^ for a detailed dis- 
cussion on polynomial generation and reproduction of subdivision schemes. 

In addition, from a practical point of view, the results in this paper, 
if combined with those in jjjj], provide a simple algebraic tool for the con- 
struction of new interesting multivariate subdivision schemes with enhanced 
properties. Note that, in the univariate case, affine combinations of existing 



schemes have been already used in [5H8|, lllHl3j for such purpose. 
The remainder of this paper is organized as follows. In section [2} we briefly 
summarize key notions concerning subdivision schemes with integer, expand- 
ing dilation matrix M. In section [31 we prove that for a non-singular multi- 
variate scalar subdivision scheme with finitely supported mask a = {a a , a £ 
Z s } and symbol a(z) = & a z a , the polynomial reproduction of degree up 

to k is equivalent to 

s ji-l 

(£Ma)(l,...,l) = |det(M)| J] Ylin-ii) and (D'a)(e) = (1.1) 

for e £ j £ Ng and |j| < k. The set S' in (II. ip is a finite set of certain 
multi-indices and r = (t±, ■ ■ ■ , r s ) £ IR 5 appears in the parametrization as- 
sociated with the subdivision scheme. The importance of condition (11.11) for 
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k = 1 is that it allows us to identify the correct parametrization of any non- 
singular or even only convergent subdivision scheme to guarantee at least 
the reproduction of linear polynomials. The parametrization determines the 
grid points to which the newly computed values are attached at each step of 
the subdivision recursion to ensure the highest degree of polynomial repro- 
duction of the scheme. Since some preliminary results discussed in jij are 
still valid for a general integer expanding matrix M, in the following we only 
prove the extensions of 0, Propositions 2.3 and 2.5]. In section HJ we use 
the results in section [3] to increase the degree of polynomial reproduction of 
existing subdivision schemes by considering their affine combinations. 

2. Background and notation 

Let M G Z sxs be an expanding dilation matrix whose spectral radius 
p(M) satisfies p(M) > 1. The r— th step of a scalar s-variate subdivision 
scheme is given by 

d£ +1) = (S a d (r) )« = Za-MpdP , d(°> G £(Z S ), r > , aGZ s , 

(2.1) 

and is defined by the finitely supported mask a = {a a e R, ct e Z s } and 
the subdivision operator 

S^.i(Z s )^£(Z s ), (5 a d) a = ^a^-Mpdp, aGZ s , (2.2) 

where £(Z S ) is the space of scalar sequences indexed by Z s . We denote the 
symbol of a scalar s-variate subdivision scheme with mask a by 

a(2) =J2 a <*z a , z = (*!, . . . , Zs ) G (C \ {0}) s , z tt = z^zT • • ■ • • C- 
Let m := |det(M)| and 

E = {e ,...,e m _i} (2.3) 

be the set of representatives of Z S /MZ S containing = (0,0,..., 0). The 
m submasks of a and the associated symbols a e (z), e G E, are defined by 
{aM Q+e , ct G Z s } and a e (^) = aMa+e 2 Ma+e , respectively. 
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Definition 2.1 (Condition Z&). We say that the mask symbol a(z) satisfies 
the zero condition of order k (Condition Z^), if 

a(l) = m and (LPa)(e) = for eeE':=E\{l} and |j| < k , 

(2.4) 

where 5 = {eq, . . . , e m _i} is i/ie sei of representatives of r L s jM Tr L s which 
contains 1 = (1, 1, . . . , 1). 

For t£| s and a 6 Z s define 

t«:=ti, r) + M-a, t^^t^-M-V, = 0, r > 0. (2.5) 

We call the sequence {t^ r \ r > 0}, t^ r ^ = {ta\ a G Z s }, the sequence of 
parameter values associated with the subdivision scheme. 

Definition 2.2. If the sequence of continuous functions {F^ r \ r > 0} with 
F ( - r \ta' > ) = da\ a e Z s , converges for all initial data d^ = {d« , ex. G 
Z s } G ^(Z/ 5 ), then we denote its limit by S^°S ^ = lim and say that S a 

is convergent. 

Definition 2.3. A subdivision scheme is called non-singular, if it is conver- 
gent, and S%°S°} = if and only if da = for all a G If . 

We denote by Hfc the space of multivariate polynomials of total degree 
k > 0. 

Definition 2.4. A convergent subdivision scheme S a with parametrization 
{t( r \ r > 0} is reproducing polynomials up to degree if for any it G H 7r 
and dW = {7r(tL 0) ), cx G Z s } the limit of the subdivision process satisfies 
S^S® = IT. 

3. Algebraic conditions for polynomial reproduction 

In this section we extend the results of Proposition 2.3 and Proposition 
2.5 in 0] to the case of a general expanding dilation matrix M G Z sxs . 

Proposition 3.1. Let S a be a non-singular s-variate scalar subdivision scheme 
that generates linear polynomials, i.e. its symbol satisfies Condition Z 2 . Then 
S a reproduces linear polynomials if and only if its parameter values are given 
by ( TO]) with 

r = mr l (D £l a(l),...,D £s a(l)) , 
where D 3 is the'yth directional derivative and is the l-th unit vector ofW. 
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Proof: By 0, Proposition 1.7] it suffices to prove the claim for the step- 
wise polynomial reproduction. Moreover, any convergent subdivision scheme 

(r) 

reproduces constants, hence we only consider the starting sequences d a = 
7r(ta' ) ) = (ta?)j, a € Z s and j = l,...,s. Then for any a £ H and 
e = (e\, . . . , e s ) £ E we get 

i(r+l) ,(r) ,(r) 

/3ez s /3gz s 



^a A//3+e (tJ, r) + M- r (a-/3)^ 

3ez s 

]T a M/ 3+e (t|, r) + M-(' +1 '(Ma + e)) - ^ a M/3+ e (M"^ +1 )(M/3 + e)) . 



(3£Z S 3 /3GZ 



= (ti, r) + M-( r+1 »(Ma + e)) - ( ]T a M/3+e (M/3 + e 

= (ti r) + M-^ +1 )(Ma + e)) - (M^ +1 )( J D £l a e (l), . . . , ^^(1)))^. 
= (ti, r) + M-( r+1 )(Ma + e) -m- 1 M-( r+1 )( J D ei a(l),...,D es a(l))) . 

where we use the fact that ajw/3+ e = a e(l) = 1 and the last equality is 
due to 0, Proposition 2.1, part (u) ] for j = €j. Thus, d^9 e is equal to 
T(t£S«) = + (M-^(Ma + e)) . = t%+ (M^Ma - r + e)) . , 
a £ Z s , if and only if Tj = m _1 D 6j a(l) for all j = 1, • • • , s. ■ 



As in [2[ we easily get the following consequence of Proposition 13.11 for 
convergent subdivision schemes. 

Corollary 3.2. Let S a be a convergent s-variate scalar subdivision scheme 
that generates linear polynomials, i.e. its symbol satisfies Condition Z 2 - Then 
S a reproduces linear polynomials if its parameter values are given by $2.5\) 
with 

T = m~ 1 (D £l a(l),...,D £ °a(l)). 
The next Proposition is crucial for the proof of our main result. 
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Proposition 3.3. Let € N, r 6 I s and q- } given by 

s ji-l 

?oO) := 1, gjOi,...,z s ) := Yl Y[(zi-£i), j = (ji,..-,J s ), zGl s . 

(3.1) 

j4 subdivision symbol a(z) satisfies 

(LVa)(l)=mq i (r), (D*a)(e) = Q e EE', j G Nq, |j| < A:, (3.2) 
z/ and only if 

{M- r py = (M^ r+1 \(x - t)Y , a G If, j G N*, |j| < k, r > . 

(3.3) 

Proof: Note that due to 0, Proposition 2.1- part (Hi) } conditions in (13.21) 
are equivalent to 

<?j( T ) = ^ a ~ Mp^-Mf,, j G Nq, |j| < k, aeZ s . (3.4) 

The proof is by induction on k. For k = we get, for any r G R s , 9o(t) — 
& a -Mf3 = (M~( r+1 ^ (at — t)) = 1. Assume next that the claim is true 

for all j 6 Nq with |j| < k — 1 and prove it for j G Nq with |j| = k. The 
polynomial gj in x of (total) degree | j | = k is of the form 

(3.5) 

£GNg, \l\<k 

Therefore, using the induction assumption and by (13.41) and (13. 5p we have 
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/3ez s pez s \£\<k 



E 


Cj >a> £ ^2 a <*-Mf3 (M r flY + 


E 


^2 a "-M/9 (M r /3) 


eew , \£\ 




feNg, |^|<fc-i 


/3GZ S 


E 


c- s , a ,i ^2 &a ~ M P (M~ r /3) f + 


E 


c w (M^+^a-r))* 


£eNg, \£\ 




£GNg, |^|<fc-l 




E 




- (M"^ +1 )( 


a - r))^ +gj(r) . 


£eNg, \e\ 


=fc V/3GZ" 







The last equality is due to the fact that 

a (r) = (q (o - M^M^+^a - r)) = ^ (M^ +1 )(a - r))" . 

£eN§, |^|<fc 

Note that all rows and columns of M r+1 , r > 0, are non-zero vectors, due 
to p(M) > 1 and det(M r+1 ) = det(M) r+1 . Thus, if all diagonal elements of 
M r+1 are non-zero, c- i a j ^ 0. If one of the diagonal elements of M r+1 is zero, 
then there exists ansxs permutation matrix P, such that for any j G Nq, 
|j| = k, the coefficient Cj )Ctj pj ^ in 

W ( ]T a a _ M/ 3 (M" r /3) £ - {M^ +1 \cx - r))M = 0. 

£ S Ng, |£|=fc \/3GZ a / 

Hence, comparing the coefficients on both sides of the above identity, we note 
that it is satisfied for all j G Nq, |j| = k, if and only if 

J2 a »^w (M- r /3) j -(M-^ +1 )(a - r)) j = 0, for all j G Ng, |j| = k. ■ 

The following result follows directly from 0, Theorem 2.6] and Proposition 
13.31 by replacing ml by Me Z sxs whenever appropriate. 
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Theorem 3.4. Let k G No- A non- singular s-variate scalar subdivision 
scheme with symbol a(z) and associated parametrization in \2. 5\) with some 
r 6K S reproduces polynomials of degree up to k if and only if 

(£> j a)(l) = m gjCr) and (D i a){e) = for e e S', |j| < k. 

We conclude by observing that, if a scheme is only convergent, then Theo- 
rem [33] leads to the following sufficient conditions for determining the degree 
of polynomial reproduction of the scheme. 

Corollary 3.5. Let k e N . A convergent s-variate scalar subdivision 
scheme with symbol a(z) and associated parametrization in A2.5\) with some 
r£R s reproduces polynomials of degree up to k if 

(£> j a)(l) = m gj(r) and (D i a)(e) = for e e E', |j|<ife. 

4. Applications and examples 

4-1. Box splines 

The results of [i[ imply that any bivariate convergent subdivision scheme 
whose symbol satisfies Condition Z 3 is an affine combination of the 3— directional 
box spline symbols from the set 

{Bj,j,h, Bj jh j, Bhjj ■ h = 0, 1, j = 3 — h} 

with 

One easily checks that e.g. the affine combination 

a(z) = 5 ■ .6221(2;) - ^212(2) - B m (z) - 2 ■ B330O) 

satisfies all conditions of Corollary 13.51 for k = 3 although none of the sum- 
mands of this affine combination does separately. The subdivision scheme 
associated with a(z) is C 1 as its second difference operator is contractive. 
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4-2. subdivision 

This example shows how to determine the degree of polynomial repro- 
duction of the approximating a/3— subdivision schemes given in [15|, page 
21] from the corresponding mask symbol a(z) instead of the iterated sym- 
bol a(z±z 2 2 , z\ z^ 1 ) ■ a(z), as it is done in [2J. We also show how to use 
affine combinations of these schemes to improve their degree of polynomial 
reproduction. 

The dilation matrix in this case is 

M =(_2 _l)> M 2 = -M, 
and, e.g., the mask symbol is given by 

a(z) = - (ziz 2 + zT 1 ^ 1 + z7 x zi + z{ 2 z 2 + ziz^ 2 + z\z 2 x ) 


+ g {zi 1 + z 2 + z 1 z 2 ~ 1 ) + - (z 2 l + z 1 + z^ 1 z 2 ) . 

The associated subdivision scheme satisfies zero conditions of at most order 2, 
see 15|. The result of Corollary 13.21 yields r = (0, 0), which implies that the 



corresponding scheme reproduces linear polynomials, if we set r = (0, 0) in 
( 12.51) . Since, the mask symbol satisfies at most zero conditions of order 2, the 
associated refmable function has approximation order 2, see [l6j]. Note that, 
similarly, the corresponding r is (0, 0) for all approximating subdivision 
schemes given in 15|, page 21]. 

Take next the symbols a,j(z), j = 1,2,3,4 of the four approximating 
subdivision schemes in [l5[ page 21] that satisfy zero conditions of order 3 
and consider their affine combination 

4 4 
j=l 3=1 

The resulting scheme still satisfies Condition Z3 and in addition the rest of 
the conditions in Corollary 13.51 for k — 3, if 

1 3 
Ai = -A 3 + A4 + 2, A2 = — -A3 — 2A4 — 1, A3, A4 G K.. 

Surprisingly, for any such choice of the parameters Ai, . . . , A 4 the resulting 
scheme is given by 

a(z) = l-^(zf 2 + ^f 2 Z2+^2+ 2 l+ 2 1^2" 2 + ^2" 2 ) + ^r 1 + ^r 1 ^2+^2 + ^l + ^l^ 2 " 1 +^ 
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and defines the interpolatory scheme considered in 15|, page 17]. Thus, by 
161 ]. this scheme has approximation order 3. 



4-3. General expanding dilation 



We use the results of [14J to define a convergent subdivision scheme with 
dilation matrix 

' 2 1 



M ={0 2)- ("> 
Similarly, we could also treat the case of any general integer expanding dila- 



tion matrix. As shown in [14j, one particular choice for the coset representa- 
tives of any such M is obtained by determining the intersection of M ■ [0, l) 2 
with Z 2 . Thus, for M in flHD we get E = {(0, 0), (1, 0), (1, 1), (1, 2)}. Ap- 
plying the strategy described in (H. llo| we get via convolution for the mask 
symbol 

eGE 

a convergent scheme associated with the symbol |a 2 (z). The symbol a 2 (z) 
satisfies ([H]) with k = 1 and E = {(1, 1), (1, -1), (-1, i), (-1, -i)}, i = 
The conditions in f!2.4[) are not satisfied for k = 2, e.g. for j = (0, 2) and 
e = (1,-1). Using the result of Corollary I3.2jw e get r = (2,1), which 
together with the zero conditions of order 2, by (l6j, implies that the scheme 
has approximation order 2. 
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